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Abstract 

In this paper the equations governing the deformations of infinitesi- 
mal (incremental) disturbances superimposed on finite static deforma- 
tion fields involving magnetic and elastic interactions are presented. 
The coupling between the equations of mechanical equilibrium and 
Maxwell's equations complicates the incremental formulation and par- 
ticular attention is therefore paid to the derivation of the incremental 
equations, of the tensors of magnetoelastic moduli and of the incre- 
mental boundary conditions at a magnetoelastic/vacuum interface. 

The problem of surface stability for a solid half-space under plane 
strain with a magnetic field normal to its surface is used to illustrate 
the general results. The analysis involved leads to the simultaneous 
resolution of a bicubic and vanishing of a 7 x 7 determinant. In order 
to provide specific demonstration of the effect of the magnetic field, 
the material model is specialized to that of a "magnetoelastic Mooney- 
Rivlin solid" . Depending on the magnitudes of the magnetic field and 
the magnetoelastic coupling parameters, this shows that the half-space 
may become either more stable or less stable than in the absence of a 
magnetic field. 

1 Introduction 

One of the main reasons for industrial interest in rubber-like materials re- 
sides in their ability to dampen vibrations and to absorb shocks. This paper 
is concerned with an extension of the nonlinear elasticity theory adopted for 
describing the properties of these materials to incorporate nonlinear mag- 
netoelastic effects so as to embrace a class of solids referred to as magneto- 
sensitive (MS) elastomers. These "smart" elastomers typically consist of 
an elastomeric matrix (rubber, silicon, for example) with a distribution of 
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ferrous particles (with a diameter of the order of 1-5 micrometers) within 
their bulk. They are sensitive to magnetic fields in that they can deform 
significantly under the action of magnetic fields alone without mechanical 
loading, a phenomenon known as magnetostriction. As a result, their me- 
chanical damping abilities can be controlled by applying suitable magnetic 
fields. This coupling between elasticity and magnetism was probably first 
observed by Joule in 1847 when he noticed that a sample of iron changed its 
length when magnetized. 

In general, the physical properties of magnetoelastic materials depend on 
factors such as the choice of magnetizable particles, their volume fraction 
within the bulk, the choice of the matrix material, the chemical processes of 
curing, etc.; see pQ for details, and also [2] for an experimental study on a 
magneto-sensitive elastomer. 

The coupling between magnetism and nonlinear elasticity has generated 
much interest over the last 50 years or so, as illustrated by the works of 
Truesdell and Toupin [3], Brown [1], Yu and Tang [5], Maugin [6], Eringen 
and Maugin [7], Kovetz |H], and others. The corresponding engineering ap- 
plications are more recent (see Jolly et al. [9j, or Dapino [lOj, for instance) 
and have generated renewed impetus in theoretical modelling (see, for exam- 
ple, Dorfmann and Brigadnov [Tlj; Dorfmann and Ogden [12]); Kankanala 
and Triantafyllidis [13j. Here, we derive the (linearized) equations governing 
incremental effects in a magnetoelastic solid subject to finite deformation in 
the presence of a magnetic field. These equations are then used to examine 
the problem of surface stability of a homogeneously pre-strained half-space 
subject to a magnetic field normal to its (plane) boundary. Related works on 
this subject include the studies of McCarthy [11], van de Yen [l5j, Boulanger 
[inillZ], Maugin [TH], Carroll and McCarthy [H] and Das et al. 

We adopt the formulation of Dorfmann and Ogden [12] as the starting 
point for the derivation of the incremental equations. This involves a total 
stress tensor and a modified strain energy function or total energy function, 
which enable the constitutive law for the stress to be written in a form very 
similar to that in standard nonlinear elasticity theory. The coupled governing 
equations then have a simple structure. We summarize these equations in 
Section [21 For incompressible isotropic magnetoelastic materials the energy 
density is a function of five invariants, which we denote here by /i and J2, 
the first two principal invariants of the Cauchy-Green deformation tensors, 
and /4, Is, Jg, three invariants involving a Cauchy-Green tensor and the 
magnetic induction vector. This formulation is similar in structure to that 
associated with transversely isotropic elastic solids (see Spencer [21])- The 
general incremental equations of nonlinear magnetoelasticity are then derived 
in Section [3l Therein we define the various magnetoelastic 'moduli' tensors 
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and provide general incremental boundary conditions. Care is needed in 
deriving the boundary equations since the Lagrangian fields in the solid and 
the Eulerian fields in the vacuum must be reconciled. 

Section m provides a brief summary of the basic equations associated with 
the pure homogeneous plane strain of a half-space of magnetoelastic mate- 
rial with a magnetic field normal to its boundary. In Section El the general 
incremental equations are applied to the analysis of surface stability. Not sur- 
prisingly, the resulting bifurcation criterion is a complicated equation, even 
when the pre-stress corresponds to plane strain and the magnetic induction 
vector is aligned with a principal direction of strain, as is the case here. The 
bifurcation equation comes from the vanishing of the determinant of a 7 x 7 
matrix, which must be solved simultaneously with a bicubic equation. To 
present a tractable example, we therefore focus on a "Mooney-Rivlin magne- 
toelastic solid" for which the total energy function is linear in the invariants 
-^1, I21 Ia, and J5. Of course, these invariants are nonlinear in the deforma- 
tion and the theory remains highly nonlinear. The bicubic then factorizes 
and a complete analytical resolution follows. In addition to the two elastic 
Mooney-Rivlin parameters (material constants), the material model involves 
two magnetoelastic coupling parameters. The stability behaviour of the half- 
space depends crucially on the values of these coupling parameters and also 
on the magnitude of the magnetic field. In particular, a judicious choice of 
parameters can stabilize the half-space relative to the situation in the ab- 
sence of a magnetic field. Equally, the half-space can become de-stabilized 
for different choices of the parameters. Thus, even this very simple model 
illustrates the possible complicated nature of the magnetoelastic coupling in 
the nonlinear regime. 

2 The equations of nonlinear magnetoelastic- 
ity 

In this section the equations for nonlinear magnetoelastic deformations, as 
developed by Dorfmann and Ogden [121 [221 [231 [21] , are summarized for sub- 
sequent use in the derivation of the incremental equations. 

We consider a magnetoelastic body in an undeformed configuration Bq, 
with boundary dBo- A material point within the body in that configuration 
is identified by its position vector X. By the combined action of applied 
mechanical loads and magnetic fields, the material is then deformed from Bq 
to the configuration B, with boundary dB, so that the particle located at X 
in Bq now occupies the position x = x{^) the deformed configuration B. 
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The function x describes the static deformation of the body and is a one- 
to-one, orientation-preserving mapping with suitable regularity properties. 
The deformation gradient tensor F relative to Bq is defined by i*" = Gradx, 
Fia = dxi/dXa, Grad being the gradient operator in Bq. The magnetic field 
vector in B is denoted H, the associated magnetic induction vector by B 
and the magnetization vector by M. 

To avoid a conflict of standard notations, the Cauchy-Green tensors are 
represented here by lower case characters; thus, the left and right Cauchy- 
Green tensors are b = FF^ and c = F^F, respectively, where * denotes the 
transpose. The Jacobian of the deformation gradient is J = det F, and the 
usual convention J > is adopted. 



2.1 Mechanical equilibrium 

Conservation of the mass for the material is here expressed as 

Jp = po, (2.1) 

where po and p are the mass densities in the configurations Bq and B, respec- 
tively. For an incompressible material, J = 1 is enforced so that p = pq. 

The equilibrium equation in the absence of mechanical body forces, is 
given in Eulerian form by 

divT = 0, (2.2) 

where r is the total Cauchy stress tensor, which is symmetric, and div is the 
divergence operator in B. The total nominal stress tensor T is then defined 
by 

T = JF'V, (2.3) 
so that the Lagrangian counterpart of the equilibrium equation (12. 2p is 



DivT = 0, (2.4) 

Div being the divergence operator in Bq. 

Let N denote the unit outward normal vector to dBo and n the cor- 
responding unit normal to dB. These are related by Nanson's formula 
nda = JF~*NdA, where dA and da are the associated area elements. The 
traction on the area element in dB may be written rnda or as T^NdA. A 
traction boundary condition might therefore be expressed in the form 

T'N = t,, (2.5) 

where ta is the applied traction per unit reference area. If this is independent 
of the deformation then the traction is said to be a dead load. 
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2.2 Magnetic balance laws 



In the Eulerian description, Maxwell's equations in the absence of time de- 
pendence, free charges and free currents reduce to 

divB = 0, cmlH = 0, (2.6) 

which hold both inside and outside a magnetic material, where curl relates 
to B. Thus, B and H can be regarded as fundamental field variables. A 
third vector field, the magnetization, when required, can be defined by the 
standard relation 

B = iio{H + M). (2.7) 

We shall not need to make explicit use of the magnetization in this paper. 

Associated with the equations (12. 6p are the boundary continuity condi- 
tions 

(B-B*) n = 0, (H -W)xn = 0, (2.8) 

wherein B and H are the fields in the material and B* and H* the corre- 
sponding fields exterior to the material, but in each case evaluated on the 
boundary dB. 

Lagrangian counterparts of B and H, denoted Bi and Hi, respectively, 
are defined by 

Bi = JF^B, Hi = F^H, (2.9) 
and in terms of these quantities equations (12. 6p become 

DivBi = 0, CnilHi = 0, (2.10) 

where Curl is the curl operator in Bq. We note in passing that a Lagrangian 
counterpart of M may also be defined, one possibility being Mi = F^M. 

The boundary conditions (12.81) can also be expressed in Lagrangian form, 
namely 

{Bi - JF-^B") ■ N = 0, {Hi- F^W) x N = 0, (2.11) 
evaluated on the boundary dBo. 

2.3 Constitutive equations 

There are many possible ways to formulate constitutive laws for magnetoelas- 
tic materials based on different choices of the independent magnetic variable 
and the form of energy function. For present purposes it is convenient to 
use a formulation involving a 'total energy function', or 'modified free energy 
function', which is denoted here by Q, following Dorfmann and Ogden [12] . 
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This is defined per unit reference volume and is a function of F and Bi: 
fl{F, Bi). This leads to the very simple expressions 

for a magnetoelastic material without internal mechanical constraints, and 

T=^-pF-\ H, = ^ (2.13) 
OF ^ ' ^ dBi ^ ' 

for an incompressible material, where p is a Lagrange multiplier associated 
with the constraint det F = 1. Note that the expression for Hi is unchanged 
except that now det = 1 in ^2. 

The Eulerian counterparts of the above equations are 

T = J-'F— H = F-'^ (2 14) 

^ dF' dBi ^ ^ 

for an unconstrained material, where = (i^^^)*, and 

where I is the identity tensor. We emphasize that the first equation in each 
of fi2.i2p - fi2.i5p has exactly the same form as for a purely elastic material in 
the absence of a magnetic field. 



2.4 Isotropic magnetoelastic materials 

In general the mechanical properties of magnetoelastic elastomers have fea- 
tures that are similar to those of transversely isotropic materials. During 
the curing process a preferred direction is 'frozen in' to the material if the 
curing is done in the presence of a magnetic field, which aligns the magnetic 
particles. If cured without a magnetic field then the distribution of particles 
is essentially random and the resulting magnetoelastic response is isotropic. 
We focus on the latter case here for simplicity, but the corresponding analysis 
for the more general case follows the same pattern, albeit more complicated 
algebraically. A general constitutive theory for the former situation has been 
developed by Bustamante and Ogden [25j and applied to some simple prob- 
lems. For isotropic materials, the energy function VL depends only on c and 
Bi® Bi, through the six invariants 

h = tr c, h = \ [(tr cf - (tr c^)] , Jg = det c = 

h = BiBi, h = {cBi)-Bu h = {c^Bi)-Bi. (2.16) 
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For incompressible materials, Is = 1 and only the five invariants Ji, I2, 
/4, J5, and Iq remain. The total stress tensor r is then expressed as 

r = -pI+2nib+2n2{hb-b'^) + 2n5B(^B + 2ne{B^bB + bB^B), (2.17) 

where Qi = dQ/dIi, and the total nominal stress tensor T as 

T = -pP-^ + 2fiiF* + 2VL2{hF^ - F%) 

+ 2n^Bi ®B + 2nQ{Bi 0bB + F^B ® B). (2.18) 

Finally, the magnetic field vector H is found from (12.151) 9 as 

H = 2{nib-^B + n5B + n6bB), (2.19) 

and its Lagrangian counterpart is 

Hi = 2{n4Bi + n5cBi + nec''Bi). (2.20) 

2.5 Outside the material 

In vacuum, there is no magnetization and the standard relation (12. 7p reduces 
to 

B* = fioH\ (2.21) 

where the star is again used to denote a quantity exterior to the material. 
Also, the stress tensor r is now the Maxwell stress r*, given by 

= fiQ^lB* ^ B* - liB"" ■ B*)I], (2.22) 

which, since divB* = and curlS* = 0, satisfies divr* = 0. 

3 Incremental equations 

3.1 Increments within the material 

Suppose now that both the magnetic field and, within the material, the 
deformation undergo incremental changes (which are denoted by superposed 
dots). Let F and Bi be the increments in the independent variables F and 
Bi. It follows from (I2.12p that the increment T in T and the increment Hi 
in Hi are given in the form 

t = AF + TBi, Hi = TF + KBi, (3.1) 
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where A,, F and K are, respectively, fourth-, third- and second-order tensors, 
with components defined by 



(3.2) 

We refer to these tensors as magnetoelastic moduli tensors. We note the 
symmetries 

— ^I3jai, ^afS — K,pa, (3-3) 

and observe that T has no such indicial symmetry. The products in (13. ip are 
defined so that, in component form, we have 

Tai = AaipjFjp + Tai/sBi^, Hi^ = T piaFifS + /Ccj/3i?;^. (3.4) 

For an unconstrained isotropic material, is a function of the six in- 
variants Ji, I2, 13, 14, 15, Iq, and the expressions (13.21) can be expanded in the 
forms 

•Aaif3j - 2^ 2^ "^""dF^dF^ ^ dF- dF„' 

m=l, m^4 n=l, 717^4 ■'^ n=l,n7i:4 ■'^ 



m=4n=l,n^4 7i=5 'Z' 

6 6 dl dl ^ d'^I 

^"/^ = E E 9# + E ^-557^' ^^-^^ 

m=4 ri=4 '/^ n=4 '/3 

where fi„ = dfl/din, ^mn = d^fl/dlmdln- Expressions for the first and 
second derivatives of J„, n = 1, . . . , 6, are given in the Appendix. 

For an incompressible material, T is given by (12. 131) i and its increment 
is then 

f = AF + TBi-pF-^ +pF-^FF-\ (3.6) 

which replaces (13.11) in this case. On the other hand. Hi is still given by 
(I2.I2P 9 and its increment is unaffected by the constraint of incompressibility, 
except, of course, since fl is now independent of Is = 1, the summations in 
equations (13. 5p omit m = 3 and n = 3. 

It is now a simple matter to obtain the incremental forms of the (La- 
grangian) governing equations. We have 

DivT = 0, DivS; = 0, Curlif^ = 0. (3.7) 
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These equations can be transformed into their Eulerian counterparts (indi- 
cated by a zero subscript) by means of the transformations 

To = J-^Ff, B^o = J-^FBu Hi, = F 'Hi (3.8) 

(with J = 1 for an incompressible material), leading to 

divf = 0, div^io = 0, curlij^o = 0. (3.9) 

Now let u denote the incremental displacement vector x — X. Then, 
F = Gradn = (grad'u)-F, where grad is the gradient operator with respect to 
X. We use the notation d for the displacement gradient gradxt, in components 
dij = dui/dxj. From fl3.8p and (13.11) we then have 

fo = Aod + ToBio, Hio = rod + KoBio, (3.10) 

where, in index notation, the tensors Aq, Fq, and Kq are defined by 

(3.11) 

for an unconstrained material. For an incompressible material J = 1 in the 
above and fl3.10p is replaced by 

f^ = A^d + T^Bi^+pd-pi, Hio = Tod + KoBio, (3.12) 

and the incremental incompressibility condition is 

divu = 0. (3.13) 

Notice that and /Co inherit the symmetries of A- and K, respectively, 
so that 

Aojisk = Aoskji, ^oij = ^oji- (3-14) 

Finally, using the incremental form of the rotational balance condition 
FT = {FTy, we find that Fq has the symmetry 

^Qijk = Fojifc, (3.15) 

and we uncover the connections 

Aojisk Aoijsk "^js^ik Tis^jk (3.16) 

between the components of the tensors Aq and r for an unconstrained mate- 
rial (see, for example, Ogden [29j for the specialization of these in the purely 
elastic case), and 

Aojisk - Aoijsk = (Tjs + p6js)6ik - (ris + p5is)6jk (3-17) 
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for incompressible materials (see Chadwick [30j for the elastic specialization). 
Following Prikazchikov [21], we decompose the tensor ^ into the sum 

Ao = A'^+J^'^ + J^'l (3.18) 

The first term does not involve any derivatives with respect to I4, Is, and 
Iq. Clearly, this term is very similar to the tensor of elastic moduli associated 
with isotropic elasticity in the absence of magnetic fields. In component form 
it is given by 

+2dikbjs^i + 2{2bijbks + SikAfjs — bjkhs — bikbjs)^2 

+2j'^{2Sij5ks - 5is5jk)^3, (3.19) 

where 

■Mj = bkkbij — bikbkj (3.20) 
and bij are the components of b. 

The terms ^0^^^^ and ^0^^^^ may be expressed in the forms 

6 

. (5) _ .0(5) ^ ^m(5) ^ 

m=l, mj^A 
6 

m=l, TTly^i 

where ^JJ-J^ = and 

-Kfiik = '^J'^cijO'sSik, -Kfiik = -^J'^iakasbij + aiajbks), 

■^Ifiik = 4:J'^{akasMj + aiajAfks), ^oj-2fc = ^Jio-kas^ij + aiajSks), 

"^Ojisk ~ OjjOjiOjsO'ki '^Ojisk ~ (fljO'i^ks ~\~ O'kO's^ij) i 

(3.22) 

with 

Tiij = aja^bik + aittkbjk, <Xi = Fi^Bi^. (3.23) 
Similarly, ^ojilfc = ^^^1 

•^ofisk ~ 2</ ^{SikTi-js + CliClsbjk + O.jO'kbis + CljCtsbik + O^iO'kbjs), 

•^ofisk ~ 4J(7Yfcs5jj + Ti-ijSks), -^ofisk ~ ^{o-iajHks + CLkds'Hij), 
Alfl, = U-'n^.-Kks. (3.24) 
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The tensor Fq is decomposed as 

To = r(^) + r'i^ + r'i^ + rf + t'^^\ (3.25) 

with components given by 

^ojik = 4:ajaiM5k + 2{aj6ik + ai6jk)^l5, 

^ojik = 4^HijM6k + 2{6ikasbjs + aibjk + djkashs + ajbik)fl6, (3.26) 

where 

■Mik = F~f^ Bi^riii + afcf2j5 + ajhjk^iQ. (3.27) 
Finally, we represent K,q in the form 

K^ = Kt'^ + K^^^ +Kf\ (3.28) 

with components 

<tj = 2JF-^{2Bi^M,, + F-^)n,, 
/eg = 2J{2aiM^j+5ijn^), 

lCfy= 2J{2akhkM^,+h,^^^). (3.29) 

For an incompressible material, the above expressions are unaltered ex- 
cept that J = 1 and all the terms ^3 and f2„3, n = 1, . . . , 6 in Aoijku ^oijk, 
fCoij are omitted. 

3.2 Outside the material 

The standard relation B = fioH in vacuum is incremented to 

B* = fioH*, (3.30) 

where B and H are the increments of B* and H*, respectively. These 
fields satisfy Maxwell's equations 

divS* = 0, curlij* = 0. (3.31) 

Finally, we increment the Maxwell stress of (12.221) to 

T* = 12^^ [B^ 0B* + B*®B* - {B* ■ B*)I] , (3.32) 

noting that divr* = 0. 
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3.3 Incremental boundary conditions 

At the boundary of the material, in addition to any apphed traction ta (de- 
fined per unit reference area), there will in general be a contribution from 
the Maxwell stress exterior to the material. This is a traction T*n per unit 
current area and can be 'pulled back' to the reference configuration to give 
a traction Jt*F~^N per unit reference area, in which case the boundary 
condition 02.51) is modified to 

T'N = Jt*F''N + ta. (3.33) 

On taking the increment of this equation, we obtain 

t^N = Jt*F-^N - Jt*F-^F^F-^ + Jt^F-^N + ta, (3.34) 

and hence, on updating this from the reference configuration to the current 
configuration, 

tin = T*n - T*d^n + {divu)T*n + 4- (3.35) 

Proceeding in a similar fashion for the other fields, we increment the 
magnetic boundary conditions (12.111) to give, again after updating, 

{Bio + dB* - {divu)B* - B*) ■ n = (3.36) 

and 

{Hio - dtW -H*)xn = 0. (3.37) 

4 Pure homogeneous deformation of a half- 
space 

Here we summarize the basic equations for the pure homogeneous deforma- 
tion of a half-space in the presence of a magnetic field normal to its boundary 
prior to considering a superimposed incremental deformation in Section 5. 

4.1 The deformed half-space 

Let Xi, X2, be rectangular Cartesian coordinates in the undeformed 
half-space Bq and take X2 = to be the boundary dBo, with the material 
occupying the domain X2 > 0. In order to minimize the number of param- 
eters, we consider the material to be incompressible and subject to a plane 
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strain in the {Xi,X2) plane. With respect to the Cartesian axes, the defor- 
mation is then defined by Xi = XXi, X2 = X^^X2, = X^. The components 
of the deformation gradient tensor F and the right Cauchy-Green tensor c 
are written F and c, respectively, and are given by 



(4.1) 



where A is the principal stretch in the Xi direction. The invariants Ji and I2 
are therefore 

Ji = /2 = 1 + A^ + X~\ (4.2) 

We take the magnetic induction vector B to be in the X2 direction and 
to be independent of xi and x^. It then follows from divB = that its 
component B2 is constant. Thus, 

Bi = 0, B2^0, Bs = 0. (4.3) 

The associated Lagrangian field Bi = F~^^B then has components 

5,1 = 0, Bi2 = XB2, 5,3 = 0, (4.4) 

and the invariants involving the magnetic field are 



Bf2, 



X-'L 



x-^h. 



(4.5) 



We may now compute the stress field using f l2.17p . (14. II) and (14.41) . The 
resulting non-zero components of r are 



Tu = 2QiX^ + 2Q2{X^ + I) - P, 

T22 = 2niX''^ + 2^2(1 + X-^) - p + 2n5X~^h + 4fi6A"^/4, 
r33 = 2r]i + 2fi2(A' + A-2)-p. 

The magnetic field H has components given by (I2.19p as 

Hi = 0, H2 = 2(^4 + X-^Q^ + X^^Q6)XBi2, = 0. 



(4.6) 



(4.7) 



Since B12 and A are constant, all the fields are uniform and the equilibrium 
equations and Maxwell's equations are satisfied. 

In view of (14. 2 p and (14. 5p . there are only two independent variables, A 
and 74. We thus introduce a specialization a;(A, J4) of the total energy Q, by 
the definition 



cu(A, h) = n{l + A" + A"", 1 + A' + A-", h, X-'h, X-^h 



(4.8) 
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from which it follows that 

uj^ = n^ + x-^n^ + A-^^e, (4.9) 

where lox = duj/dX,uJi = du/dl4,. Hence, 

Til - T22 = Xux, H2 = 2XBi2UJi. (4.10) 

4.2 Outside the material 

From the boundary conditions (12.81) applied at the interface X2 = X2 = 0, 
we have B2 = B2 and = = 0, while from fl2.2ip it follows that 
B* = B^ = and H2 = yUo ^-^2 = l^o^B2- Outside the material we take 
the magnetic field to be uniform and equal to its interface value, Maxwell's 
equations are then satisfied identically, B* therefore has components 

B{ = {], Bl = B2 = X-^Bi2, B* = 0, (4.11) 

and H* has components 

Ht = 0, H^ = fio'B2=fio'^''Bi2, H; = 0. (4.12) 

From these expressions, we deduce that the non-zero components of the 
Maxwell stress f l2.22p are given by 

r*u = -T22 = -WbI = -1/^0 'A-2/4 = r3^3- (4.13) 

The applied mechanical traction on X2 = required to maintain the plane 
strain deformation has a single non-zero component r22 — T22. 

5 Surface stability 

We now address the question of surface stability for the deformed half-space 
by establishing a bifurcation criterion based on the incremental static solution 
of the boundary-value problem. Biot [32] initiated this approach, which has 
since been successfully applied to a great variety of boundary- value problems; 
see Ogden [33] for pointers to the vast literature on the subject. 
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5.1 Magnetoelastic moduli 

First we note that since Fij = for i ^ j and Bn = B13 = several simpli- 
fications occur in the expressions for the components of the magnetoelastic 
moduli tensors A^, Tq, Kq. In particular, we have 

Aoiijk = 0, JCoij = 0, for J ^ k, 

ToiiS = Tosii = Foiil = Tqih = 0, 

Toijfe = 0, for 2 i^hi-i. (5.1) 
For subsequent use we compute the quantities 

O = -4.01212, 26 = »4oilll + »4o2222 — 2^01221 — 2^01122, C = »4o2121, 
A— ro211, e = ro222 — roil2, f — ^011, 9' — ^022- (5-2) 

Explicitly, we obtain 

a = 2A2(1]i + fis) + 2/41^6, 

b={X^ + X-^){Qi + Q2) + h[>^~^^5 + (6A-^ - 2)Qe] 
+ 2(A^ + A-^ - 2)(Qn + 21^12 + ^^22) 

+ AhiX-' l)[ni5 + ^25 + 2A-2(l]i6 + 1^26)] 

+ 2/2A-^(fi55 + 4A-2fi56 + 4A-^1]66), 

c = 2A-2(ni + 1^2) + 2J4[A-2(]5 + (2A-^ + 1)1^6], 

c^ = 25,2A[A-2(]5 + (A-^ + l)l]6], 

e = 4Bi2\-'[n5 + 2A-2fi6 + (1 - A^)(fii4 + ^24) 

+ (A-2 - X^){Q^, + O25) + - l)(l^i6 + ^26) 

+ 74(^254 + A '^^55 + 2A ^1^46 + 3A ^fige + 2A ^r2ge)], 

/ = 2(A"^1]4 + + A^l^e) 

y = 2(A^1]4 + + A-^l^e) + 4/4(A^l]44 + 2f]45 + 2X-^Q46 

+ A-'i755 + 2A-^(]56 + A-^1^66). (5.3) 
In terms of the energy density a; (A, I4) we have the connections 

a — c — Xcux, 2(6 + c) — }?uj\\i e — —2Bi2}?uj\a-i g — 2A^(a;4 + 2I4U44), 

(5.4) 

where uw = d'^uj/dX^, lox^ = d'^u/dXdh and a;44 = d'^u/dlj. 

5.2 Incremental fields and equations 

We seek incremental solutions depending only on the in-plane variables Xi 
and X2 such that U3 = and 5^03 = 0. Hence Ui = Ui{xi,X2) and BiQi — 
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BiQi{xi,X2) for i = 1,2 and p = p{xi,X2). In the following, a subscripted 
comma followed by an index i signifies partial differentiation with respect to 

Xi, t 1,2. 

The incremental version (13.131) of the incompressibility constraint reduces 
here to 

Mi,i + M2,2 = 0, (5.5) 
and hence there exists a function ip = iIj{xi,X2) such that 

ui = i^,2, U2 = (5.6) 
Similarly, equation (13.91) -:! reduces to 

^^01,1 + 5(02,2 = 0, (5.7) 

and the function = 0(xi,X2) is introduced such that 

Bioi = 0,2, Bio2 = -0,1- (5.8) 
The incremental equations of equilibrium (I3.9p i simplify to 

^011,1 + ^021, 2 = 0, Toi2,l + To22,2 = 0. (5-9) 

From the identities (15. ip . the only non-zero components of the incremental 
stress To are found to be 

Ton = (Ann + p)uia + An22^i2,2 + 5o2ron2 - P, 

Tq21 = (wAo2112 + P)U2,1 + ^2121^1,2 + -Bo^ro211, 
^012 = (Al221 + P)ui^2 + Al212^^2,l + -Bo^roi21, 

To22 = (A2222 + P)U2,2 + A22nMl,l + ^0310222 " P- (5.10) 

Also, equation (13.91) 9 reduces to 

Hm,2-Hio2,i=0, (5.11) 
wherein are the only non-zero components of Hi, which, from (15. ip . are given 

by 

Hm = roi2l(Ml,2 +"^2,1) + l^OuBm, H102 = roil2Ml,l + ro222W2,2 + A^022-B«02- 

(5.12) 

In terms of the functions ip and equations (15.91) and (15. lip become 

(Ann ~ A1122 — Al22l)'0,112 + A212l'0,222 ~ rO1120,ll + rO1210,22 = P,l, 
(A2222 — A1122 — Al22l)^,122 + Al212'0',lll — (r0121 — rO222)0,12 = ~P,2, 
(r0222 — r0112 — r012l)^,112 + rO121'0,222 + ^^0220,11 + AI^O110,22 = 0. (5.13) 
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We eliminate p from the first two equations by cross-differentiation and ad- 
dition and obtain finally the coupled equations 

+ 2fcV',1122 + C'?/',2222 + (e - rf)0,112 + #,222 = (5-14) 

and 

dll),222 + (e - d)%l)^ii2 + /0,22 + #,11 = (5.15) 

for ip and 0. 

5.3 Outside the material 

In vacuum, Maxwell's equations (13.311) hold for B and H. From the second 
equation, and the assumption that all fields depend only on xi and X2, we 
deduce the existence of a scalar function 0* = 0*(xi,X2) such that 

Ht = -ct>% H* = -<P% Hl = Q. (5.16) 
Equation (13.301) then gives 

B{ = -^io'P% B*2=-f^o<P% ^3 = 0, (5.17) 
and from (I3.3ip i we obtain the equation 

0:11 + 0:22 = (5.18) 

for 0*. Finally, the incremental Maxwell stress tensor (13.321) has non-zero 
components 

= X-'Bi2<j)% = f3*3 = -f*2, f*2 = -A-iS^20:i = f*i. (5.19) 

5.4 Boundary conditions 

We now specialize the general incremental boundary conditions of Section [331 
to the present deformed semi-infinite solid. First, for ta = 0, the incremental 
traction boundary conditions (I3.35P reduce to 

T02I + T^iU2,l - = 0, to22 + T^2U2,2 " = 0, (5.20) 

on X2 = 0. Putting together the results of this section, using (14.131) . (15.21) . 
O, dSISD, fCTJj) . flETUj) and dSinD, we express the two equations fCTJj) as 

(^22 + 1/^0 '^"'^4 - C)^,ll + C^,22 + rf0,2 + >^-' Bi2(l)\ = 0, (5.21) 
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and 

(26 + c-r22 + i/io'A-2/4)V',ii2 + cV',222 + e<^,ii + #,22-A-'S;2(/'*i2 = 0, (5.22) 

which apply on 0:2 = 0. In obtaining the latter we have differentiated (I5.2UI) 9 
with respect to Xi and made use of fl5.13l) i . 

Next, the incremental magnetic boundary conditions 03.36P and f l3.37p 
reduce to 

5^02 + B1^U2,2 - = 0, Hm - H*U2,i - Ht = (5.23) 

on a;2 = 0. Using again the results of the preceding sections, we write these 
as 

A-^5;2^,i2 + 0,1 - /io0:2 = 0, (5.24) 

and 

(^0 ^-^^2 - C?)^,ll + #,22 + /0,2 + 0^1 = (5.25) 

on X2 = 0. 

5.5 Resolution 

We are now in a position to solve the incremental boundary value problem. 
We seek small-amplitude solutions, localized near the interface X2 = 0. Hence 
we take solutions in the solid {x2 > 0) to be of the form 

i/j = Ae-^'''^e'^''\ = kDe-^'''^(^^''\ (5.26) 

where k > Q (2,71 /k is the wavelength of the perturbation) and s is such that 

3?(s) > (5.27) 

to ensure decay with increasing X2 (> 0). 

Substituting (15.261) into the incremental equilibrium equations (15.141) and 
(I5.15p . we obtain 

(cs^ - 2bs^ + a)A - s{ds^ + d-e)D = 0, 

s{ds^ + d-e)A-{fs^ - g)D = Q. (5.28) 

For non-trivial solutions to exist, the determinant of coefficients of A and D 
must vanish, which yields a cubic in s^, namely 

{cf -d^)s^ -[2bf + cg + 2{d-e)d]s^ + [2bg + af -{d-ef]s^ -ag = 0. (5.29) 
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From the six possible roots we select Si, S2, S3 to be the three roots satisfying 
(I5.27p . We then construct the general solution for the solid as 

3 3 
ip = J2 ^j-e-'^'^^^^e'^^S (p = kJ2 Dje-^'^''^(^^''\ (5.30) 
i=i i=i 

where Aj, Dj, j = 1,2,3, are constants. 

For the half-space 0:2 < (vacuum) we take a solution 0* to (15.181) that 
is localized near the interface X2 = 0. Specifically, we write this as 

0* = ifcCV^^e'^^^S (5.31) 

where C* is a constant. 

The constants Aj and Dj are related through either equation in (15.281) . 
From the second equation, for instance, we obtain 

Sj{ds'^ + d — e)Aj + {fs^- — g)Dj = 0, j = 1, 2, 3; no summation. (5.32) 

We also have the two traction boundary conditions (I5.2ip and (I5.22p . which 
read 

(c - - \^lQ^\'^h){Al + A2 + A3) + c{s\Ai + slA2 + ^^^3) 

- d{siDi + S2l^2 + 53/^3) - X-^BnC- = 0, (5.33) 

and 

(^22 - \^lQ^\-^h - 26 - c){siAi + S2A2 + S3A3) 

+ c{slAi + slA2 + slA^i) + (e - dsl)Di 

+ (e - dsl)D2 + (e - dsl)D^ - X-^BnC- = 0. (5.34) 

Finally, the two magnetic boundary conditions (I5.24p and (15.250 become 

\-'Bi2is,A, + S2A2 + S3A3) - (D, + D2 + Ds) + fioC" = 0, (5.35) 

and 

{d - i2o^X-^Bi2){Ai + A2 + A3) + d{slAi + slA2 + slA^) 

- /(si/^i + 52/^2 + 53/^3) -C* = 0. (5.36) 

In total, there are seven homogeneous linear equations for the seven un- 
knowns Aj, Dj, j = 1, 2, 3, and C*. The resulting determinant of coefficients 
must vanish and this equation is rather formidable to solve, particularly since 
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it must be solved in conjunction with the bicubic (15.291) . It is in principle pos- 
sible to express the determinant in terms of the sums and products S1 + S2+S3, 
S1S2 + S2S3 + S3S1, S1S2S3, and to find these from the bicubic (15.291) . similarly 
to the analysis conducted in the purely elastic case (see Destrade et al. [3l] ) . 
However, the resulting algebraic expressions rapidly become too cumbersome 
for this approach to be pursued. 
Instead, we propose either 

(a) to turn directly to a numerical treatment once Vt has been determined 
by curve fitting from experimental data for a given magnetoelastic solid, 

or 

(b) to use a simple form for that allows some progress to be made. 

Regarding approach (a), we remark that, as emphasized by Dorfmann 
and Ogden [121 ESI ES, Elj, there is a shortage of, and a pressing need for, 
suitable experimental data and for the derivation of functions Vt from such 
data. In the next section we focus primarily on the analytical approach (b). 

5.6 Example: a "Mooney-Rivlin magnetoelastic solid" 

As a prototype for the energy function fi, we propose 



n = i/i(0)[(l + 7)(/i - 3) + (1 - 7)(/2 - 3)] + fi,\ah + f^h), (5.37) 



where /i(0) is the shear modulus of the material in the absence of magnetic 
fields and a, f3, 7 are dimensionless material constants, a and [3 being mag- 
netoelastic coupling parameters. For a = /3 = 0, (I5.37P reduces to the strain 
energy of the elastic Mooney-Rivlin material, a model often used for elas- 
tomers. 

In respect of (I5.37P the stress r in (I2.17P reduces to 



T = -pI+ i/io(l + 7)& + |/^o(l - l)ihb - b") + 2/io ® B, (5.38) 



Clearly, equation (15.380 shows that the parameter a does not affect the stress. 
By contrast (3, if positive, stiffens the material in the direction of the magnetic 
field, i.e. a larger normal stress in this direction is required to achieve a given 
extension in this direction than would be the case without the magnetic 
field. On the other hand, by reference to (I5.39p . we see that a provides a 



while H in (12.190 becomes 




(5.39) 
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measure of how the magnetic properties of the material are influenced by the 
deformation (through b). If /5 = the stress is unaffected by the magnetic 
field. On the other hand, if a = then the magnetic constitutive equation 



f l5.39p is unaffected by the deformation. Thus, a two-way coupling requires 
inclusion of both constants. 

The quantities defined in fl5.2p and (15.31) now reduce to 

a = /i(0)A^ 2b = fi{0){X^ + X'^ + pX-^h), c = fi{0){X-^ + f3X-^h), 



d = ^Jfi^^fi{0)(3X-'Bl2, e = 2^ fi^' fi{0)(3X-'Bi2, 

f = fi^\aX~' + g = fi^\aX^ + (5.40) 

where B12, a dimensionless measure of the magnetic induction vector ampli- 
tude, and I4 are defined by 



B12 = 5z2/v//io/i(0), h = Bl. (5.41) 

Note the connections 

26 = a + c, e = 2d. (5.42) 
Now we find that the bicubic fl5.29p factorizes in the form 

(s^ - 1) (s^ - A^) [aA^ + pX^ - (a + pX^ + aph)s^] = 0, (5.43) 

and it follows that the relevant roots are 



^1 = 1' ^2 = A^ .3 = AW //'^t^.f - (5-44) 

a -I- pX^ + aph 

Note that for S3 to be real for all A > and all B12, the inequalities 

a>0, p>0 or a>0, P>0 (5.45) 

must hold. (The case in which there is no magnetic field corresponds to 
a = (3 = 0.) It is assumed here that these inequalities are satisfied, so that 
S3 is indeed a qualifying root satisfying (15.271) . 
The equation (15.320 becomes 

Sj{s^j - l)pX-^Bi2Aj - [(aA-2 + /5)s| - aX^ - P]Dj = 0, j = 1, 2, 3, (5.46) 

where 



Aj = i//io V(0)Ai, = fi^^Dj, (5.47) 
and the sj are given by (15.440 . 
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Next, consider the four remaining boundary conditions fl5.33p - fl5.36p . In 
order to keep the number of parameters to a minimum (so far, we have A, B12, 
a, f3), and to make a simple connection with known results for the surface 
stability of an elastic Mooney-Rivlin material, we assume that there is no 
applied mechanical traction on the boundary X2 = 0, and hence 

r22 = T^2 = K'A-2J4. (5.48) 
The boundary conditions f l5.2ip -f l5.25p now read 

[I + {(3- l)/4](ii + A2 + is) + (1 + f3h){slAi + slA2 + slA;) 

- pXBi2{siDi + S2D2 + s^Ds) - \B12C* = 0, 
(A^ + 2 + 2/3/4) (siii + S2A2 + S3A3) - (1 + f3h){slAi + slA2 + SIA3) 

+ (3XBi2[{sl - 2)z3i + {si - 2)D2 + {si - 2)D,] + XB12C' = 0, 
Bi2isiAi + S2A2 + ssAs) - \{Di + D2 + D3) + XC" = 0, 
\Bi2{P - l)(ii + A2 + is) + \PBi2islAi + SIA2 + slAs) 

- (a + P\^){s,£)i + S2D2 + ss^s) - A'C* = 0. (5.49) 

From the seven equations fl5.46p and (15.490 , we have derived a bifurcation 
criterion (vanishing of the determinant of coefficients) using a computer alge- 
bra package, but it is too long to reproduce here. It is a complicated rational 
function of the four parameters A, B12, a, (3. However, it is easy to solve 
numerically, and for the numerical examples we fix the material parameters 
a and (3 and find the critical stretch Acr in compression as a function of B12. 
For B12 = 0, we recovered the well-known critical compression stretch for sur- 
face instability of the elastic Mooney-Rivlin material in plane strain, namely 
Acr = 0.5437 [32], as expected. For Figure la (Figure lb), we set a = 0.5 
(a = 2.0) and curves for /? = 0.0, 0.5, 1.0, 1.5, 2.0 are shown. We found that 
Acr is an even function of B12 and we therefore restricted attention to positive 
B12 (within the range < B12 < 3). The behaviour as B12 becomes larger 
and larger (not shown here) indicates that the half-space becomes more and 
more unstable in compression. Moreover, it can even become unstable in 
tension (Acr > !)• The figures also clearly demonstrate that for some values 
of a, f3, and B12 the critical stretch ratio is smaller than that for the purely 
elastic case (Acr < 0.5437), in which cases the magnetic field has a stabilizing 
effect. 

Turning back to a phenomenological approach, we remark that the en- 
ergy function (15.370 has quite good curve-fitting qualities for moderate fields. 
There are four parameters at hand, namely /i(0), a, /3, 7, two of which, /i(0) 
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and (3, may be determined from shear tests. Indeed Dorfmann and Ogden 
|24j show that in general the shear modulus for isotropic nonlinear magne- 
toelasticity is 2[Qi + ^2 + -^4^5 + /4^^6(3 + 2k^)], where k is the amount of 
shear in a simple shear test. Here the modulus is independent of k, and is 
given by 

fi{Bi2) = m + 2/Xo i/?/4. (5.50) 

This highlights the role of (3 in increasing the mechanical stiffness of the 
material — through the shear modulus. Jolly et al. [9J conducted double 
lap shear tests on magneto-sensitive elastomers containing 10, 20, and 30% 
by volume of iron particles. From their Figure 7, we see that in the range 
< B12 < 0.5 Tesla, the variations of fi{Bi2) resemble those of a parabolic 
profile such as the one suggested by (15.501) . For the 10% iron by volume 
elastomer specimen. Table 1 in Jolly et al. [H] gives /i(0) = 0.26 MPa, and 
at B12 = 0.5 Tesla, we read off their Figure 7 that /i(0.5) - /i(0) ~ 0.07 MPa, 
indicating that (3 ^ 0.18. Similarly, for the 20% and the 30% iron by volume 
elastomer specimens we find f3 ^ 0.53 and (3 ~ 0.72, respectively. 

Figure 2a (Figure 2b) illustrates the variation of the critical compression 
stretch with the amplitude of the dimensional magnetic induction vector, 
from to 0.5 Tesla, for the 20% (30%) iron by volume elastomer, and for 
several values of a. We remark than the presence of the magnetic field 
makes the two specimens slightly more stable than in the purely elastic case 
because all the critical compression stretch values are smaller than 0.5437. It 
is also clear that increasing the value of a makes the half-space more stable. 
However, it is worth noting that the 30% iron by volume specimen is slightly 
less stable than the 20% iron by volume specimen for the same values of a. Figure 2 



A Derivatives of the invariants with respect 
to F and Bi 

We derive the expressions for the first derivatives of the six invariants with 
respect to F, 

2F- ^ = 2(c F -c F) 



dF dF 



dh . dh „ dL 



2^3^aO TTT^ = 0, —— = 2Bi^{Fi^B, 



dF °" dF dF- 



dh 

-- 2{FiyBiyCapBip + Fi^CypBij^Bia)., (A.l) 



dF 
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and with respect to Bi, 

=0, ^ = 0, ^ = 0, 



dBia dBia dBia 

— = 2Sia, = 2Ca/3Bip, — — ^ = 2Ca'yC^pBip. (A.2) 

OJJia Onia OtSia 

The second derivatives of the invariants are computed as follows: first, 
the second derivatives with respect to F, 



dFiadFjp 

dFiadFjp 
h 



dFiadFj/j 



= 0, 

= 26ijBiaBii3, 



dFiadFj/3 

OF aF ^ '^[^ijiCayBiyBip + CpyBiyBia) + 5apFiyBiyFjsBis 

+ FiyBiyFjaBip + FjyBiyFiijBia + bijBiaBifs]; 

(A.3) 

next, the mixed derivatives with respect to F and Bi, 



dFi^dBi^ ' 9F,«aS,, ' dFiM, ' 9F,c.9S^^ 



5 



dFi^dBi^ 



If-h 



dFi^dBi^ 



'^^apFiyBi^ + 2Bi^Fii3, 



'^FijjCayBi^ + 2FiyBi^Cai3 + 2FiyCyi3Bi^ + 26ai3FiyCysBis; (A.4) 



finally, the second derivatives with respect to Bi, 



dBiadBip ' dBiadBifj ' dBi^dBifj 

5^/4 <9^-?^5 <9^-^6 

— — — — = 25^/3, o J-, = 2cq^, — — „ = 2ca^c^^. (A. 5) 

oUiaOnifj oniaOnifj oniaO-oip 
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Figure Captions 



Figure 1. Dependence of the critical stretch Acr < 1 for instabihty in com- 
pression for a magnetoelastic Mooney-Rivhn sohd in plane strain on the 
non-dimensional measure B12 of the magnetic field for several values of the 
magnetoelastic coupling parameters a and /?. 

Figure 2. Dependence of the critical stretch Acr < 1 for instability in com- 
pression for a magnetoelastic Mooney-Rivlin solid in plane strain with the 
dimensional measure B12 of the magnetic field, for several values of the mag- 
netoelastic coupling parameters a and (3. 
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Figures 
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(a) a = 0.5 (b) a = 2.0 

Figure 1: Dependence of the critical stretch Acr < 1 for instabihty in com- 
pression for a magnetoelastic Mooney-Rivhn sohd in plane strain on the 
non-dimensional measure B12 of the magnetic field for several values of the 
magnetoelastic couphng parameters a and [3. 
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Figure 2: Dependence of the critical stretch Acr < 1 for instabihty in com- 
pression for a magnetoelastic Mooney-Rivhn sohd in plane strain with the 
dimensional measure B12 of the magnetic field, for several values of the mag- 
netoelastic coupling parameters a and j3. 
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